Gold fragmentation induced by stopped antiprotons 
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A natural gold target was irradiated with the antiproton beam from the Low Energy Antiproton 
Ring at CERN. Antiprotons of 200 MeV/c momentum were stopped in a thick target, products 
of their annihilations on Au nuclei were detected using the off-line 7-ray spectroscopy method. In 
total, yields for 114 residual nuclei were determined, providing a data set to deduce the complete 
mass and charge distribution of all products with A > 20 from a fitting procedure. The contribution 
of evaporation and fission decay modes to the total reaction cross section as well as the mean mass 
loss were estimated. The fission probability for Au absorbing antiprotons at rest was determined 
to be equal to (3.8 ± 0.5) %, in good agreement with an estimation derived using other techniques. 
The mass-charge yield distribution was compared with the results obtained for proton and pion 
induced gold fragmentation. On the average, the energy released in p annihilation is similar to that 
introduced by ~1 GeV protons. However, compared to proton bombardment products, the yield 
distribution of antiproton absorption residues in the N-Z plane is clearly distinct. The data for 
antiprotons exhibit also a substantial influence of odd-even and shell effects. 

PACS numbers: 25.43.-|-t, 25.85.Ge 



I. INTRODUCTION 



The large energy of almost 2 GeV released in nucleon- 
antinucleon annihilation had awoken hopes to observe 
some unique nuclear reactions induced in this way. Espe- 
cially energetic antiprotons were presumed to reach the 
deep interior of the nucleus. Exotic processes, like phase 
transitions to quark-gluon plasma and explosive decay of 
that hot system had been expected to occur ||, but 
were not verified in experiments performed shortly after 
the commissioning of the Low Energy Antiproton Ring 
at CERN 1^, ^. Nevertheless, the character of the re- 
actions starting with antiproton absorption in nuclei is 
quite unique in comparison with reactions induced with 
protons or heavy ions. Whereas the excitation energy 
carried in by post-annihilation mesons is quite large, the 
linear and angular momentum transfer as well as the mat- 
ter compression are reduced, particularly in the stopped 
antiproton case which is preceded by the exotic atom 
phase. Hence, one may investigate a clear thermal reac- 
tion aspects with suppressed collective motion complica- 
tions. 

Such phenomena were intensively studied during the 
LEAR era for stopped and for energetic antiprotons. The 

f ectra for neutrons and light charged particles 
|, |lO|, mass yield distributions 0, ^,jr 
and characteristics of the fission fragments [Jl 



were measured for a wide range of targets, 
excitation energy derived from these studies. 
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and w 300 MeV for heavy nuclei absorbing stopped and 
1.2 GeV antiprotons, respectively, compares well with the 
average values obtained for protons which have about 1 
GeV larger energies, i.e. approximately the nucleon rest 
mass. 

Yield distributions of residual nuclei were studied as 
soon as more intense p beams were provided by LEAR. 
Many targets were irradiated with antiprotons, mainly at 
rest energy, but only few of them were examined in detail: 
"^*Cu PI, 92^95 Mo jnl 98mo (i|, "'^'Ag H, "'^^Ba [||, 
^^^Ho and ^^-'^Ta ]l(|. Average quantities such as the 
mean mass removed from the target as well as individ- 
ual yield features, e.g. isomeric ratios, were investigated. 
Clear odd-even effects in the mass and charge yield dis- 
tribution were observed. Theoretical calculations, based 
on intranuclear cascade -|- evaporation models, were able 
to reproduce only the gross features and failed to predict 
yield dependence on the detailed N and Z . The odd- 
even phenomenon, although postulated to be present and 
used to model the yield distribution for years pO| , pT| and 
sometimes reaching very large values |2^, still seems to 
be almost completely unexplored theoretically in a more 
quantitative way. 

One of the most distinct features observed in reactions 
with stopped antiprotons is the large probability {~ 0.1) 
of very small energy transfer, when the target nucleus 
looses only one nucleon in annihilation and is left ex- 
cited below the nucleon separation energy. Nuclear spec- 
troscopy studies of the relative yield of both types of 
these residual {At — 1) nuclei (a neutron or a proton lost 
in such soft antiproton absorptions) were used to estab- 
lish a new and powerful method of probing the nuclear 
periphery composition |2^, 24 , 2^, |26| . 
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The irradiation of the heavy, gold target gave us a 
chance to study the competition between evaporation 
and fission induced by antiproton absorption. The yield 
distribution of heavy residual nuclei complements the re- 
sults obtained from on-line measurements of neutron and 
charged particle spectra Since gold is a commonly 
used target, there was also the opportunity to compare 
these data with a rich set of information gathered from ir- 
radiations with energetic protons or pions. In particular, 
the yield distribution after the reaction of 1 GeV protons 
with Au has been extensively studied in older and recent 

and also with 
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7-ray spectroscopy measurements 
a new method using the mass-charge spectrometry for in- 
verse kinematic reactions Our preliminary data 
of the reaction p + Au at rest were already published 
, this work presents the complete results obtained af- 
ter fitting the Y{A, Z) yield distribution. 

The article is organized as follows. In Section II wc 
briefly present some details of the experiment, in Sec. Ill 
data analysis is described together with the yield fitting 
method. Experimental results are presented and initially 
discussed within Sec. IV. In Sec. ^ we compare the re- 
sults of this experiment with those obtained with other 
projectiles impinging on Au and with data coming from 
studies of antiproton absorption in various targets. Fi- 



nally, our main conclusions are presented in Sec. VI 



II. EXPERIMENT 

A thick target of natural gold was irradiated with the 
antiproton beam of 200 MeV/c momentum from LEAR 
facility. The target of the total thickness of 549 mg/cm^ 
was composed of ten foils of 80, 30, 30, 30, 2, 30, 30, 
37, 80 and 200 mg/cm^, starting from the beam side. 
The initial energy of the antiprotons, equal to about 21 
MeV, was reduced in the scintillation counter (from Pi- 
lot B) and in some additional moderators (mylar, silicon) 
to about 6.5 MeV at the target surface. Such an ar- 
rangement assured that the majority of antiprotons was 
stopped in few central Au foils. The very central and 
extremely thin foil of 2 mg/cm^ was applied to monitor 
the X-ray activity, while the last and thickest one (200 
mg/cm^) was used to check the secondary reactions level. 

Two scintillation counters, SI and S2, were used to 
control beam intensity and transverse dispersion. The 
first anticounter SI had a hole of 10 mm diameter and the 
active area of counter S2 was a disc of the same diameter. 
Consequently the signal S1S2 indicated particles going 
towards the target. The irradiation lasted 15 minutes 
with the total number of antiprotons equal to (9.25 ± 
0.35) X 10^ (STS2 number). 

Monitoring of the target activity started 13 min after 
the irradiation and continued at CERN during one week; 
afterwards the spectra were collected in Warsaw, the last 
one was taken more than a year after target activation. 



Two HPGe detectors were used at CERN, a 7-ray counter 
of 15% relative efficiency for all foils and an X-ray counter 
for the thinest one. In Warsaw two more efficient 7- 
ray detectors were applied, of 20% and of 60% relative 
efficiency, and a third X-ray detector for the thin foil. 

All collected spectra were analyzed with the program 
ACTIV 32, Gamma-ray lines were identified by 

their energies, half- lives and intensity ratios. The de- 
cay data were taken from the eighth edition of the Table 
of Isotopes Q. 

Experimental yields for all detected residual nuclei, 
normalized to 1000 p stopped in the target, are listed 
in Table |. The independent yields represent the total 
number of nuclei, summed over all isomers. Cumulative 
yields include also the yields of all /3-decay precursors of 
a given isotope. Besides that, there are presented partial 
yields for some isomers not representing the whole pro- 
duction for a given (A,Z) pair as well as some production 
limits for Hg nuclei. Mercury may be produced from gold 
after p absorption in a charge exchange reactions, when 
one of the annihilation tt"*" , ir^ pions exchanges charge 
with a target neutron. Such a phenomenon was ob- 
served for some targets irradiated with stopped antipro- 
tons Ijl^, |l^ ^ [1^, where nuclei of target charge plus 
one were produced at level ranging from 0.5 to 5 per 1000 
p. On the other hand, for some other targets, studied in 
the neutron halo project, rather low upper limits (0.5- 
2%) were given |g5[. Our data obtained for Au, except 
for the i^^Hg isomers with low 7 intensity, indicate that 
such an effect should happen very rarely. 

Initially, the distribution of the activity induced in in- 
dividual target foils was estimated several hours after 
the irradiation with the use of the measurement of the 
i86ir(g) yield in each foil. Five inner foils (30, 30, 2, 30, 
30 mg/cm^) gathered about 90% of the total activity and 
only these foils were monitored later. This reduced the 
7-ray self-absorption effect. The distribution of the tar- 
get activity was determined more precisely afterwards on 
the basis of yields obtained for six evaporation residues: 

186lr, 184jj.^ 183Qg^ ISlj^g^ 157i3y ^^^^ 152j3y^ 

age, five inner foils stopped (89.7 ± 2.1)% of the whole 
number of antiprotons. These data, together with the re- 
sults for ^^^Au and -'^^^Au, were used also to estimate the 
yield introduced by secondary reactions with particles 
(mainly pions and neutrons) produced after antiproton 
annihilation on target nuclei. It was done by comparing 
the number of given nuclei produced per foil thickness 
unit, averaged for three inner foils, with the similar re- 
sult obtained for the last, thickest foil. The upper limit 
for the secondary reactions leading to ^^^Au is equal to 
3%, the limit for -'^^^Au is about 2.6% and for the rest of 
quoted isotopes it does not exceed 2%, i.e. in all cases it 
is below the yield uncertainties. The negligible influence 
of the secondary reactions on our results is additionally 
confirmed by very small upper limit given in Table ^ for 
^^^Au, a (n,7) reaction product. 
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TABLE I: Experimental and fitted yields of residual nuclei from gold fragmentation induced with stopped antiprotons. Yield 
type: I - independent, C - cumulative, in brackets - total yield fitted for isotopes for which only one isomeric state was observed. 



Nuclide 


Tl/2 


Experiment 


Fit 


Type 


Nuclide 


Half-life 


Experiment 


Fit Type 




[N/lOOOp] 








[N/lOOOp 






Au* 


2.7 d 


< 0.6 






16b rri 

im 


30.1 h 


14.4 ± 0.8 


14.2 


C 


196 A m2 

Au 


9.7 h 


1.52 ± 0.18 






163 rp 

im 


1.81 h 


10.1 ±0.6 


10.1 


G 


196 A 

Au 


6.2 d 


75.0 ± 3.4 






161 rP™ 

im 


33.0 m 


8.6 ± 1.1 


8.6 


G 


195TT„m 

Hg 


4i.6 li 


< 3.1 






161 T7„ 

Er 


3.2 h 


2.59 ± 0.57 


2.73 


T 
i 


195tt„k 
Mg'^ 


9.9 h 


< 11.7 






160 

Er 


26.6 li 


9.8 ± 0.8 


9.9 


c 


195 A 

Au 


186 d 


38.7 ±3.1 






159t? 

Er 


36.0 m 


8.2 ± 0.8 


8.2 


G 


194 A „ 

Au 


38.0 h 


11.8 ±0.6 






157t, 


8.2 h 


8.0 ± 0.4 


7.9 


G 


193TT„m 


11.8 li 


< 0.5 






155 Ti , 

Dy 


9.9 h 


6.2 ± 0.3 


6.2 


G 


193TT„g 


S.o n 


< 2.8 






i D 


C O J 

5.3 d 


0.27 ± 0.09 


0.30 


T 
i 


193 A 
Au 


17.7 h 


9.9 ± 2.2 






154tt„k 

Ho** 


11.8 m 


r> r\0 1 r\ OA 


(3.16) 


G 


192tt„ 

Hg 


4.9 h 


< 1.3 






153rpi 

i b 


2.3 d 


4.77 ± 0.24 


4.78 


G 


192 A 

Au 


4.9 h 


8.5 ±0.6 


8.37 




Gd 


242 d 


0.13 ± 0.04 


0.12 


T 
i 




73.8 d 


1.61 ± 0.07 


2.6 




152^-^ 

Uy 


2.4 h 


o oci 1 nil 

3.50 ±0.11 


3.60 


G 


19lTT„m 

Hg 


51.0 m 


< 1.7 






152iTii 

i b 


17.5 h 


r\ CI 1 ri 1 o 

0.57 ± 0.13 


0.67 


T 
i 


195tt„k 


49.0 m 


< 2.8 






151 rpl 

1 b 


17.6 h 


3.55 ± 0.20 


3.55 


G 


191 A ., 

Au 


3.2 h 


6.7 ± 0.4 


6.8 




151 1 


124 d 


A _L A AO 

0.27 ± 0.08 


0.24 


T 
i 


191 

l^t 


2.9 d 


7.9 ±0.8 


8.6 




150 T-v 

Uy 


7.2 m 


1 AC _l_ A O /I 

1.95 ± 0.34 


1.95 


G 


190 TT 

Hg 


20.0 m 


< 2.2 






150rni s 

i b*^ 


3.5 h 


2.85 ± 0.27 


(2.9) 


C 


190 A 

Au 


42.8 m 


5.0 ±0.4 






149 J 

Gd 


9.3 d 


3.12 ± 0.20 


3.13 


c 


190jj.m2 


3.3 h 


2.04 ±0.15 J 






148fTTi e 

i b* 


1.0 h 


1.13 ± 0.15 


(1.73) 


c 


190 T 

Ir 


11.8 d 


1.95 ± 0.10 






I47ri J 
Gd 


38.1 h 


2.9 ±0.2 


2.9 


G 


189 dj. 


10.9 h 


12.5 ± 1.1 


12.2 




147t:i,, 

Eu 


O /I 1 J 

24.1 d 


C\ A A _l_ A OA 

0.44 ± 0.20 


0.51 


T 
i 


189 T„ 

Ir 


13.2 d 


5.1 ±0.8 


4.8 




146^1 J 

Gd 


48.3 d 


O O /I 1 A 1 O 


2.23 


G 


188 TDi 

Ft 


10.2 d 


13.2 ± 0.4 


10.6 




146 ti,, 

Eu 


4.6 d 


A "7 1 A AO 

0.73 ± 0.08 


0.71 


G 


188t„ 

ir 


41.5 11 


5.6 ± 0.4 


5.8 




145t7,, 

Eu 


5.9 d 


o 1 c 1 n o /I 

2.15 ± 0.34 


2.19 


G 


Ft 


2.4 h 


9.8 ± 1.5 


10.2 




143 ID™ 

Fm 


265 d 


2.0 ± 0.4 


1.9 


G 


Ir 


10.5 h 


7.3 ± 1.0 


6.8 




139 ri„ 

Ge 


138 d 


1.2 ±0.1 


1.2 


G 


18613+ 
Ft 


2.0 h 


10.4 ± 1.8 


• fi 




135 ri„ 
Ge 


17.7 h 


A 1 1 A 1 

U.oi dz U.zi 


0.78 


G 


186T„m 

Ir 


2.0 li 


< 4.0 






132 T „ 

La 


4.8 h 


A OO 1 A AO 

0.32 ± 0.08 


0.34 


T 
i 


ir"^ 


16.6 h 


7.8 ± 0.8 






131t3„ 

Ba 


11.5 d 


C\ CO 1 A O /I 

0.53 ± 0.24 


0.62 


G 


185t„ 

Ir 


14.4 h 


16.3 ± 1.1 


16.5 


U 


129 ri„ 
Gs 


32.1 h 


A aCv _L A OA 

O.oO ± 0.20 


0.56 


T 
i 


186 rv„ 
Us 


93.6 d 


5.2 ±0.8 


4.1 


T 
1 


Ae 


36.4 a 


A C "1 _l_ A AC 

0.51 ± 0.05 


0.52 


G 


184 Tjj. 

Ft 


17.3 m 


8.2 ± 0.8 


9.0 




124t 
i 


4.2 d 


C\ "X O 1 A AT 

0.15 ± 0.07 


0.15 


T 
i 


184t„ 

Ir 


3.1 n 


8.3 ± 0.5 


9.6 




121rp„ 

ie 


16.8 d 


A O A 1 All 

0.39 ±0.11 


0.38 


G 


183t^ 

Ir 


58.0 m 


14.2 ± 2.3 \ 




u 


103 T3„ 

Hu 


39.3 d 


A /I O 1 A CV7 

0.48 ± 0.07 


0.48 


G 


183 /-^„m 

Us 


9.9 h 


11.4 ±0.8 J 




Kj 


ic 


4.3 d 


A O O _L A A CT 

0.23 ± 0.05 


0.23 


T 
i 


I83r\„g 
Us* 


13.0 h 


9.7 ±0.6 




u 


JN b 


35.0 d 


A _l_ A AT 

0.d5 ± 0.07 


0.65 


T 
i 


183 D„ 

He 


70.0 d 


2.22 ± 0.54 


2.23 


T 
1 


95 rv^ 


64.0 d 


Air* 1 A AO 

0.15 ± 0.03 


0.16 


G 


182t„ 

Ir 


15.0 m 


11.8 ± 0.6 


14.8 




93 TV T„m 

Mo 


6.9 h 


A O O 1 A A /I 

0.23 ± 0.04 


(0.23) 


T 
i 


I82n„ 
Us 


22.1 h 


9.8 ± 1.0 


9.9 


T 
i 


At 


3.3 d 


A r O 1 A A C 

0.58 ± 0.05 


0.59 


G 


181d„ 

He 


19.9 h 


24.9 ± 1.4 


25.5 


u 




83.4 d 


A OA _L C\ C\ A 

0.30 ± 0.04 


0.30 


G 


180 rv„ 
Us 


21.7 m 


19.9 ±0.8 


19.9 


u 


88v 
Y 


1 IT? J 

107 d 


A O "1 _l_ A "1 "1 

0.81 ± O.ii 


0.83 


T 
i 


He 


19.5 m 


20.8 ± 2.0 


20.3 




87v 
Y 


3.3 d 


A 1 A A /i 

0.75 ± 0.04 


0.74 


c 


He 


13.2 m 


15.5 ± 2.4 


15.7 


c 




14.7 h 


A O O 1 A A /I 

0.38 ± 0.04 


0.39 


G 


w 


2.3 h 


20.8 ± 1.2 


21.2 




85 Q 

or 


64.8 d 


A O 1 1 All 

U.oi zb U.li 


0.83 


G 


ia 


56.6 h 


3.4 ± 1.1 


2.1 


T 
i 


84 -nl- 

Hb 


32.8 d 


A TO _L A 1 O 

0.78 ± 0.12 


0.81 


T 
i 


176m„ 

ia 


8.1 h 


26.5 ± 2.0 


26.4 


u 


83 Til 

Hb 


86.2 d 


A OA _l_ C\ "X C 

0.89 ± O.io 


0.85 


G 


175 rp„ 

ia 


10.5 h 


21.2 ± 1.5 


21.6 




82 T) urn 

Hb 


6.5 h 


A /t A 1 All 

U.4U zb U.ii 


(0.43) 


T 
i 


175 -rrr 

Hi 


70.0 d 


2.3 ±0.6 


1.6 


I 


Hr 


35.3 h 


A O 1 1 A AO 

0.31 ± 0.08 


0.30 


I 


173rn 

Ta 


3.1 h 


16.1 ±2.7 


15.5 


c 


75 n 

Se 


112 d 


0.36 ± 0.05 


0.36 


c 


173Hf 


23.6 h 


5.8 ±0.7 


5.7 


I 


^^As 


17.8 d 


0.51 ±0.08 


0.51 


I 


i^^Lu 


1.37 y 


0.54 ±0.13 


1.3 


I 


^2 As 


26.0 h 


0.22 ± 0.07 


0.22 


I 




36.8 m 


8.1 ±0.7 


8.1 


c 




14.1 h 


0.27 ±0.07 


0.27 


I 


172 jjf 


1.9 y 


12.5 ± 1.4 


12.3 


I 


69 


13.8 h 


0.34 ±0.04 


(0.32) 


I 


"^Lu 


6.7 d 


1.34 ± 0.20 


1.36 


I 


5«Fe 


44.5 d 


0.36 ±0.07 


0.38 


c 


i^^Lu 


8.24 d 


20.8 ±1.0 


20.9 


c 


48y 


16.0 d 


0.19 ±0.03 


0.19 


c 


170jjf 


16.0 h 


16.1 ±1.4 


16.1 


c 


^«Sc 


83.8 d 


0.24 ± 0.07 


0.24 


I 
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TABLE I: - continued. 



Nuclide 


Half-life 


Experiment 


Fit 


Type 


Nuclide 


Half-life 


Experiment Fit 


Type 






[N/lOOOp] 










[N/lOOOp] 




^'^"Lu 


2.0 d 


3.8 ±0.7 


3.9 


I 


^^Ar 


1.8 h 


0.22 ±0.05 0.22 


C 


le^Lu 


34.1 h 


16.0 ±0.6 


16.0 


C 


^^Na 


15.0 h 


0.27 ±0.07 0.27 


C 


169Yb 


32.0 d 


2.15 ±0.67 


2.18 


I 


^^Na 


2.6 y 


<2.7 0.31 


C 


166 Yb 


56.7 h 


15.5 ±0.8 


15.5 


C 


■^Be 


53.3 d 


< 3.5 — 


C 



III. DATA ANALYSIS 

The method based on ofF-hne 7-ray spectroscopy has 
some limitations. The most important is the necessity 
to use some phenomenological model to reconstruct the 
yields of unobservable products ^ 0. In 

similar experiments with other projectiles the data are 
sparsely spread over the N-Z plane. In this case less de- 
tailed models may be applied and the results obtained 
are only a first order approximation of the true yield dis- 
tribution Y{A,Z), even when cumulative cross sections 
are involved in the fitting procedure. Moreover, the pre- 
cision of tabulated absolute 7-ray intensities sometimes 
leaves much to be desired and, finally, the decomposi- 
tion of spectra with a few hundreds of lines, as they are 
measured for heavier targets, becomes a challenge for the 
persistence of the evaluator. 

Targets irradiated with antiproton beams, much less 
intense than the proton beams, have a rather low activ- 
ity level. In particular, yields obtained for the fission 
fragments were close to our detection limit. Keeping in 
mind one of our goals, the estimation of the probability 
for antiproton induced fission, the data analysis needed 
special care and some feedback. A primary set of yield 
results obtained from the spectrum analysis served as an 
input for the model distribution fitting at its early stage, 
when the best approach was searched for. This relates 
to the choice of the final formula as well as to the divi- 
sion of the data to subsets assuring the lowest total x^- 
Afterwards, a modeling procedure was applied to check, 
confirm or eliminate some doubtful experimental yields. 
For some mass regions it appeared necessary to apply an 
additional or separate evaluation, and we describe it at 
the end of this section. 



A. Fitting procedure 

The formula, used to describe the yield distribution 
was rather complex in order to be as universal as possible 
and to test various models. This complexity mainly arose 
from the aim of taking into account cumulative yields 
and from introducing the {N, Z) evenness corrections. 
The general formula was factorized into two components, 
mass and charge distributions, Ya and Yzp, respectively. 



YiA,Z)=YAxYzp. 



(1) 



The distribution over the mass (the main distribution 
ridge) was modeled with the exponential of a fourth-order 



polynomial with parameters ai — a^, 



(2) 



This was useful for testing the fit in broader mass regions, 
where the ridge shape may change more rapidly. 

The form of the second factor in Eq. (|^), Yzp, (the 
charge distribution) was multiplied by the odd-even cor- 
rections Fp. When needed, an additional component, 
containing the sum of yields of the decay predecessors of 
given {A, Z) isotope, was added here 



Yzp ^2^Fp{Z + ck,N-ck) x e 

k=0 



-(Z+cfe-Zp)™/2(T^ 



(3) 



where the term with fc = corresponds to the indepen- 
dent Yzp{A, Z) yield and the terms with fc = 1 -^5 stand 
for the precursors contributions. The upper limit of the 
sum over k was set to 5, because the charge distribution 
for given A is rather narrow and neglecting fc > 5 did not 
change the sum by more than 1%. 

The most probable charge path Zp and the charge dis- 
tribution width a were expressed as a third-order poly- 
nomials of A, with parameters ag — ag and aio — 013, 
respectively, 



- arA + asA'^ + agA^, 



a — aio + aiiA + 012^^ ± aisA^ 



(4) 



(5) 



The value of the factor c in the sum of the yield cumu- 
lation for a given A depended on the side of the stability 
valley on which the given isotope lies, 

{1, EC, decay (neutron-deficient nuclei); 
— 1, /3~ decay (neutron- rich nuclei). 

(6) 

The power index w in the exponent argument in Eq. 
(^ was allowed to be different for Z > Zp and Z < Zp, 



ai6i Z > Zp 
an, Z < Zp, 



(7) 



where aig was always set to 2 and 017 = 2 or 017 = 1.5 was 
used in order to test the asymmetric charge distribution 
in the latter case. 
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FIG. 1: The path of the most probable atomic number 
corresponding to the best fit values for the parameters as — ag 
in seven mass regions (solid lines). Positions of nuclei, for 
which a cumulative or an independent yield was determined, 
are denoted by circles and triangles, respectively. Dashed 
line - valley of stability approximated by the relation Z/j = 
A/(1.98 + 0.0155A^/^) ||]. Open square indicates the target 
nucleus, vertical and horizontal dotted crosses are plotted for 
the magic N — 82 and the closed shell Z = Q4 and for the 
magic TV = 50 and the closed shell Z = 40. 



Finally, the odd-even correction was assumed to be a 
simple factor depending on the Z and TV evenness com- 
bination 



Fp - 



Z - even, N - even; 
Z - odd, N - even; 
Z - even, N - odd; 
ai4ai5, Z - odd, N - odd. 



(8) 



A division of the whole A-Z plane into subregions may 
be treated as yet another model parameter. To have 
control over it we have plotted positions of all detected 
nuclei in the N-Z plane; this appeared to be very helpful 
for a preliminary determination of the Zp path, espe- 
cially for regions with many data. The final mass region 
division is illustrated in Fig. |l] by the solid lines show- 
ing Zp fitted for seven data regions. Mass range limits 
were fixed to get the smallest total for the whole data 
range and to possibly simplify the model for the course 
of the Zp path. We have tested many alternative divi- 
sions, especially for the region of the heavy evaporation 
residues (143 < A < 183). The fitting applied to broader 
A ranges than those listed in Table |l| resulted in at least 
one order of magnitude larger values, mainly due to 
rapid changes in the Zp course at A ^ 162 and A — 150. 
For three separate regions of the lighter products with A 
= 121-139, 82-103 and 24-75 the limits were defined by 
the grouping of the experimental data. 

As may be seen from Table |l[ presenting the final pa- 
rameters, sometimes the best fit is obtained when the 
number of parameters exceeds the number of data. This 
was done by fixing some parameters when the others were 



fitted, and vice versa. Various combinations and order of 
fixing/releasing of parameters as well as their total num- 
ber were tested. The shape of the mass ridge could be 
parametrized with a maximum of four parameters, the 
most probable Z path was approximated via a parabola 
except for two cases and the charge distribution width 
was constant or changed linearly with mass. Odd-even 
corrections were applied only for three heaviest mass re- 
gions, where a larger number of points and smaller rel- 
ative errors of the experimental data allowed to get a 
reliable fit. The shape of the charge distributions was 
modeled better with the use of an asymmetric form for 
the evaporation residues, lying further from the stability 
valley (see Fig. ^. For fighter, fission products the Zp 
path goes closely along the valley of stability and here a 
symmetric Gaussian shape was more adequate. 

The normalized charge distributions for six mass inter- 
vals are plotted in Fig. || against the normalized charge 
difference {Z — Zp)/a. This reduces the distributions to 
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FIG. 2: Normalized Z distributions for six mass ranges. The 
fitted function (line) is the exponent exp{ — \{Z — Zp)/a\™), 
where w is the charge distribution power index, Zp the most 
probable Z and a the Z distribution width. Experimental 
values (A - independent, ■ - cumulative) are normalized as 
described in the text. All distributions have the same normal- 
ized width equal to unity (please note the asymmetric shapes 
for four distributions of heavier isotopes, with different slope 
for both sides). 
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TABLE II: Best fit parameters obtained for tiie yield distribution model applied to six mass regions. 



Paramater Mass range 





163-182 


150-161 


143-149 


121-139 


82-103 


24-75 


Ya ai 


-93.98(2) 


-90.63(3) 


-479.2(1) 


18.13(9) 


-30.97(5) 


-2.07(11) 


CL2 


1 0951 f 11 


1 060 (2) 


6 486(1) 


-0 207(11 


690fll 


-0 0059^181 




-0.00310(1) 


-0.00302(2) 


-0.02192(3) 


-0.000706(5) 


-0.00390(6) 


0.00029(3) 


(24 






-0.3(2)E-7 


0.884(4)E-5 






7i n,- 


9R 9fi('4-~l 


1 36 fifl 

-LOU .L* 1 i- J 




-7 9('Vi 


-142.1(1) 


0.77(10) 


(27 


-0 fifiSI (9\ 


-1.201(1) 




595('3l 


4 S86('1 \ 


456('2~l 


(28 


0.003081(1) 


0.00485(5) 


-0.000329(4) 


-0.00029(2) 


-0.0456(1) 


-0.000329(2) 


(2g 








0.8(14)E-7 


0.000153(1) 




a aio 


-0.82(2) 


1.17(4) 


1.12(8) 


0.91(9) 


1.15(4) 


2.12(11) 


an 


0.0105(2) 










-0.013(2) 


Fp ai4 


0.67(3) 


0.82(5) 


0.58(7) 


1. 


1. 


1. 


ai5 


0.74(3) 


0.91(6) 


0.90(8) 


1. 


1. 


1. 


w ai6 


2. 


2. 


2. 


2. 


2. 


2. 


an 


1.5 


1.5 


1.5 


1.5 


2. 


2. 


X^/NDF 


0.045 


0.006 


0.164 


0.131 


0.082 


0.020 



the same width in the case where a is not constant in the 
given region. The fitted function is the simple exponent 
eyi'p{—\{Z — Zp) / ) , then, for comparison, experimental 
data Y{A, Z)e are normalized with three factors coming 
from the fit 

Y{AZ)e 

^ Y{A)xFp{N,Z)xfj- ^' 



Here, the mass distribution Y{A) is as in Eq. (§),the odd- 
even correction Fp {N, Z) follows Eq. (||) and the factor 
// corresponds to the independent yield fraction in the 
case of the cumulative yield 



r.Z+0.5 

Fp{Z,N) I e'^^-^p^^/^a^dz 

Z-0.5 



[Fp{Z + ck,N -ck) 



k=0 



Z+ck+a.5 

X / e 

Z+ck-0.5 



(10) 



-(z+cfe-Zp)"'/2<T^ 



dz] 



r 



Sometimes charge distributions are normalized to unity 
integral over Z to get the total yield for a given A equal 
directly to Ya |3^, p9|. However, when odd-even correc- 
tions are used the distribution of the total yield cannot be 
described by a simple continuous function as in Eq. (H). 
Also the generalization of the charge distribution shape 
with the two-valued (or released) index w leads to prob- 
lems in obtaining an analytical form of the normalization 
factor for this function. Hence our distributions for six 
mass regions are normalized only to the same width a = 
1, not to the same integral. 



B. Treatment of the heaviest residues 

It is a well known fact that none simple phenomenolog- 
ical model can properly describe the charge distribution 
of yields for isotopes lying near the target nuclide ||2^, EgI . 



The main reason for this is the asymmetric, non-Gaussian 
shape of the charge distribution, with the most probable 
Z and width rapidly changing with mass. Such a phe- 
nomenon is clearly seen when one uses a longer section of 
the Y{A, Z) yield, along constant Z value instead of con- 
stant A. Fig. H presents isotopic Ya{A) distributions ob- 
tained in this experiment for elements ranging from Au to 
Lu. Also isotopic distributions obtained for the heaviest 
element after stopped antiproton absorption on ^^^Yb, 
^^^Nd and ^^°Te targets 35 exhibit such a behavior: a 



steep and narrow distribution for the target element Zt, 
a flat and broad distribution for the Z^-l element and 
deformed quasi-Gaussian shapes for some smaller Z ^ with 
the deformation on the heavy mass side decreasing with 
increasing distance from Zt- Even though the low mass 
side for all elements but Zt may be described with the 
same slope, the slope at the higher mass side changes 
rapidly and cannot be fitted well with a fixed isotopic 
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distribution asymmetry, i.e. with the unique, constant 
ail parameter. As a consequence, the heaviest elements 
should be excluded from the global fit and their Ya{A) 
yields have to be fitted separately for a given Z. 

170 IBO 190 170 180 190 170 180 190 
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FIG. 3: Isotopic yield distributions for the nine heaviest 
residual elements produced in the fragmentation of gold with 
stopped antiprotons. Dashed lines show yields obtained from 
the fit applied to the region f 63 < A < f 82; triangles, full and 
open circles represent the experimental yields: independent, 
cumulative and cumulative transformed to independent, re- 
spectively. Solid lines illustrate the corrections of the isotopic 
yield distribution for the heaviest isotopes of the elements 
from Pt to Hf, and for the lightest Au isotopes (see Sec. ^ 
for discussion). 

The method of the yield completion for the heaviest 
elements is recursive: at the beginning we estimate the 
lacking yields of the lighter Au isotopes. With the use 
of these results cumulative, experimental yields for Pt 
are converted to independent ones and the isotopic dis- 
tribution for this element is evaluated. Then, a similar 
procedure is applied for Ir, Os and so on. The method 
was applied down to Ta and Hf elements, where the 
yields for A > 175 were corrected. Finally, the summed 
Y(A) yields for A > 175 presented with a line in Fig. ^ 
are the combination of results of both evaluations: the 
163 < A < 182 region global fit and the procedure de- 
scribed above. 

The platinum distribution was the most laborious case. 
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FIG. 4: Mass distribution of residual nuclei produced by the 
reaction of stopped p with Au. The cumulative sum of ex- 
perimentally observed yields for a given mass is shown with 
crosses, the corresponding fitted Yi^A) yield for all input data 
(averaged over Z for given A) is represented by full circles. 
The line shows the sum of fitted yields completed with inter- 
polated yields for mass intervals not fitted. For the heaviest 
isotopes {A > 175) an additional correction of the mass yield 
data was performed, as described in Sec. Ill , 



due to lack of radioactive isotopes above mass 191 and 
owing to the strong odd-even N effect (up to 30%), ob- 
served for this even-Z element produced after low-energy 
absorption of antiprotons. The overall shape of the iso- 
topic distribution was assumed to be similar to that ob- 
served for Tm residues after Yb fragmentation with an- 
tiprotons p5[ , with an increasing enhancement of yields 
for three heaviest nuclei, which for Pt are these with the 
mass numbers 194, 195 and 196. The number for ^^^Pt 
obtained in this way (16/lOOOp) was confronted with the 
result of another estimation, based on the so-called halo 
factor dependence on separation energy Sn of the neu- 
tron from target nucleus p5[ . For heavy nuclei with Sn 
close to 8 MeV the halo factor is of the order of 4 5, 
hence, using the ^^^Au yield the ^(^^^Pt) should be be- 
tween 14 and 23 per lOOOp. These two estimations are 
consistent. For the lightest Pt isotopes we assumed that 
the steep slope coming from the 163 < A < 182 region fit 
is a good approximation, this assumption was used also 
for consecutive, lower Z elements (we have checked the 
justification of this approach using the lightest mass data 
points of these isotopes in an additional test fit). At last, 
on the basis of observed changes of Pt yield for odd and 
even N isotopes an appropriate correction was applied. 



C. Region of the a-decay 

There is a narrow bump at A = 147 in the mass yield 
distribution, a feature observed for Au ||2^, |2^ and for 
Ta M, Eo| target fragmentation after reactions with pro- 
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tons, heavy ions and stopped antiprotons. The enhance- 
ment of the cross section in this region was suggested 
to be the result of a-decay of nuclei above the = 82 
shell but the authors abstained to estimate quanti- 
tatively this effect due to lack of charge dispersion curves 
which could not be fitted for limited experimental data. 
We have done such an estimation for our data perform- 
ing a preliminary fit for isotopes not affected by a-decay, 
i.e. for 154 < A < 161, adding ^^"Dy with an inde- 
pendent yield. Taking into account a charge dispersion 
yield Y{A, Z) obtained in this way we have calculated 
appropriate decay corrections for experimentally mea- 
sured yields of ^^^Tb (-^2.9% correction), ^^^Gd (+2.8%), 
i52Dy (+7.2%), 152Tb (+6.0%), ^^iTb (+14.8%), ^^^Gd 
(+13.6%), 150Tb (+32.5%), i49Gd (+5.5%), ^^^Gd (- 
13.7%), i47eu (-12.1%) and i^egd (-23.7%). The ex- 
perimental results listed in Table || are the corrected ones 
and were used for fitting in two mass ranges affected by 
this effect. As can be seen from Fig. ||the corrections ob- 
tained are too small to remove the local yield maximum 
at A = 147 (crosses show the yield before correction, cir- 
cles after that). Therefore, a-decay alone cannot explain 
fully such a feature and the observed yield enhancement 
in this region should be partially ascribed to the closed 
TV = 82 shell influence. 



IV. RESULTS 

The experimental data are presented in Table |, to- 
gether with the fit results for isotopes representing full 
yield for a given pair of A and Z. Results are normalized 
to yield per 1000 p with the total number of antiprotons 
stopped in the target (9.25 x 10*). The final mass yield 
distribution is presented in Fig. The cumulative sum 
of all yields observed for a given mass number A is here 
compared with the total yield Y{A) obtained from the fit 
via summation of all fitted Y{A, Z) values over Z, or from 
the interpolation between fitted mass regions. The global 
curve of the fitted Y{A) yield, when compared with the 
summed experimental yields, forms its exact skyline in 
almost the whole region of the evaporation residues. A 
deviation from this rule is observed for three mass ranges: 
the heaviest, with A > 176, a few mass numbers around 
A = 147 and all fission fragments {A < 120). Except 
for the second region (affected by the a-decay), this is 
the result of prevailing accumulation of the isobaric yield 
by non-detectable isotopes. The depression of observed 
yield of the heaviest evaporation residues is narrow but 
deep, with a maximum decrease to about 40% of the fit- 
ted Y{A) for A w 192, and comes mainly from the stable 
Pt isotopes produced. For fission products, where the 
Zp path goes over the stability valley, the observed yield 
is strongly suppressed and its outline reaches only about 
20-50% of the fitted yield. 

Leaving out two heaviest masses, the maximal Y{A) 
yield is reached at mass 180 but the largest individual 
Y{A, Z) production is fitted for i^^W. The small yield 



peaks observed for some even masses {A = 180, 176, 170, 
...) are due to the strongest odd-even effect for some even 
Z lying almost on the Zp path. On the other hand, the 
global mass yield minimum appears between A = 105 
and A — 120. As numerous 7-lines of strongly populated 
heavier nuclei covered this region, no valuable production 
limits can be given for this region and we have to stick 
to the interpolated curve. 

After evaluation of the mass yield curve it is possi- 
ble to estimate the relative yields for different reaction 
channels. The fission fragments mass range should be 
treated with some care as their multiplicity is equal to 
2 or greater when one takes into account any multi- 
fragmentation process. Assuming that all residues with 
AO < A < 120 are binary fission products (i.e. two 
heavy residues per antiproton) and neglecting the light- 
est A < AO masses, we have obtained the summed fission 
yield. Comparing this number with the total Y{A) inte- 
gral in the mass limits from 40 to 196 we have extracted 
the probability of gold fission induced with stopped an- 
tiprotons to be (3.8 ± 0.5)%. The lighter mass region 
(yl = 10 — 40), not taken into account in fission due to 
possible multiplicity > 2 and/or the not fully negligible 
chance to have a fission partner in the A > 120 region, 
constitutes additionally less than 0.9% of the total yield 
(the error quoted above takes this into account). Our re- 
sult compares well with the fission probability of 3.1(3)% 
obtained in an experiment where fission fragments yields 
were measured with PIN diodes |4^ and is substantially 
larger then the value of 1.5% derived from another ex- 
periment using also on-line technique 143]. 



V. DISCUSSION 
A. Antiprotons versus other projectiles 

The properties of reactions induced by stopped an- 
tiproton absorption can be investigated by comparison 
with yield distributions obtained for other, more "clas- 
sical" projectiles. We have confined this comparison 
to the gold target as the literature is quite rich here 

^ti, 123, r 
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|, IJ, The two other popular 

neighbor-mass targets, Pb and Ta, represent rather dif- 
ferent decay scenarios, with respectively more and less 
pronounced fission channel. 



1. Mass yield curve 

First, we present a rather qualitative comparison with 
the yield curve shapes extracted for protons. Figure 
shows the summed isobaric Y{A) yield obtained for 
stopped antiprotons plotted together with yield distri- 
butions resulting from Au fragmentation by 0.49, 0.8, 
1.0 and 3.0 GeV protons H, |o). Since the yields 
for stopped antiprotons and protons are measured in dif- 
ferent units, we have normalized our yield axis with an 
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FIG. 5: Comparison of the mass yield distribution obtained 
for protons and stopped antiprotons. Curves for protons at 
0.49 GeV, 1.0 GeV and 3.0 GeV adopted from pj, for the 
inverse kinematic reaction of Au on H at 800 MeVfrom |2E| , 
|3o|| . The yield axis for antiprotons was normalized with a 
factor of 1.75. 



FIG. 6: Average mass removal from the Au targets irradiated 
with: pious; •, protons (Refs. cited in Table |l|) and 
stopped antiprotons (pion data are shifted right by the tt rest 
mass). The hatched band shown for antiproton data reflects 
the error in AA and a possible range of energy deposition for 
stopped p annihilation with one of the target nucleons. 



arbitrary factor equal to 0.57, providing the concordance 
between p's and 1 GeV protons results in the 150-170 
mass range. It should be stated here that the yield curve 
presented for fission residues in the case of 1 GeV protons 
was fitted with only 5 mass points [pTf . 

The most striking differences between the curves shown 
in Fig. 1^ are seen for the fission region. The fission prob- 
ability for gold excited by protons, estimated as in Sec. 
i^, is equal to « 6.5%, « 3.7% and 3.3% for 1 GeV, 0.8 
GeV and 0.49 GeV protons, respectively. Then, 800 MeV 
protons seem to correspond to stopped antiprotons, but 
fission takes only a small part of the total yield and the 
comparison of distribution shapes in the evaporation re- 
gion is much more adquate. Such inspection leads to the 
conclusion that stopped antiprotons match with protons 
at 1 GeV. 

Besides the level and width of the fission hump, a sec- 
ond major feature distinguishing the mass yield shapes 
observed for stopped antiprotons and protons is the dis- 
tribution for the heaviest masses close to the target. 
Here the experimental situation is much better than in 
the fission case: more reliable nuclear spectroscopy data 
is additionally confirmed by the results of inverse kine- 
matic measurements. The yield distribution for protons 
is rather unchanged in energy range from 0.5 to 1 GeV 
and forms a plateau between A w 175 and A « 194. On 
the contrary, for antiprotons in this mass range, not only 
for the gold target ||l5|, ^ |3^ , the yield slowly decreases 
from A « At— 20 to A oi (At— 3) and then strongly rises, 
reaching an absolute maximum at A = (At — 1). The 
enhancement of yield for few masses closest to the target 
may be explained by two mechanisms. The first one is 
the soft antiproton absorption, where almost all annihi- 
lation pions miss the rest of the target nucleus. Then the 



(At — 1) nuclei are left with very low excitation energy. 
Only after antiproton absorption on nucleons occupying 
a deeper states rearrangement of nucleon configura- 
tions results in a mass loss of one or two additional units. 
The probability of the production of {At — 1) nuclei is 
quite large, about 10% for targets used in our nuclear 
periphery studies and the results obtained for 

gold are also of this order of magnitude (cf. Table |). 
The second mechanism leading to low excitation ener- 
gies is the class of all processes where the annihilation 
pions escape unabsorbed by the target nucleus but due 
to the sizeable total 7r-nucleus cross section excite this 
nucleus enough to emit few nucleons. The probability 
of such kind of quasi-elastic meson escapes may be quite 
large when the nuclear diffuseness and partial opacity are 
considered [ p7[ . 

A quantitative comparison between antiprotons and 
other projectiles is presented in Table III and Fig. ||. 
For this purpose we have calculated the average mass 
removed from the target, AA, defined as: 



AA = A*^- 



Y{A)AdA 



(11) 



Y{A)dA 



where equals At for protons and pions or At — 1 for 
antiprotons. Amm is the lower integration limit adjusted 
to get all single heavy residues and Amax is equal to At 
for protons and pions or to At — 2 for antiprotons. The 
residual mass At — 1 is ignored in the integration of the 
reaction yields for antiprotons since it attests no reaction 
(soft p absorption). 

To get consistent and comparable results, for each data 
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set taken from the literature we have appHed a uniform 
method to determine AA. Only experimental data repre- 
senting the highest (approximately the whole) cumulative 
yield for a given A were used to construct the mass yield 
distribution curve. The absolute errors of the quantities 
presented in Table III were estimated to be of the order 
of 10-20%, but the relative errors should be smaller. In 
addition to proton data, results for pions absorbed by the 
gold target are presented, their energy range is limited 
as compared to the rest of data but coincides with the 
kinetic energy of pions emitted in antiproton annihila- 
tion. The average mass removal from the target nucleus 
smoothly correlates with the projectile energy. As can 
be seen from Table III and Fig. ^, the average mass 
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such an inspection for some cases quoted in Table III 
bring differences into prominence, we have recalculated 
Zp{A) to its distance 6 stab from the line modeling the 
beta stability valley (defined as in the caption to Fig. 
|l|). Such a presentation was earlier applied to study the 



Au + p, 0.8 GeV 
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p + Au, 11.5 & 300 GeV 

p + Au, 3 & 6 GeV 

p + Au, 1 GeV 
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removed from the gold target by stopped antiprotons is 
only slightly lower than AA obtained for 1 GeV protons. 



2. Charge distribution 

The information on the reaction mechanism, obtained 
from the investigation of mass yield distributions, may 
be enriched by the examination of other features of the 
Y{A, Z) yield topography. Two such properties were 
compared for results of gold fragmentation by antipro- 
tons and other projectiles: the course of the fitted Zp 
path and the charge dispersion width. Figure ^ illustrates at s 




FIG. 7: Upper part: Distance of the fitted most probable Z 
from the beta stability line: protons/heavy ions compilation 
at saturation limit |^] , antiprotons at rest energy compilation 
Au at 800 MeV on H 



protons at 3/6 GeV 



39|, |§, protons at 11.5/300 GeV 
and fission for protons at 1 GeV 
Lower part: Charge distribution widths as a function of 
residue mass. The results of the same evaluations as in the 
upper part are presented (with unchanged notation), except 
for protons at 1 GeV and at 3/6 GeV, for which the width 
parameter is equal to that obtained for 11.5/300 GeV [^ . 



TABLE III: Characteristics of the mass yield distribution af- 
ter the reaction of different projectiles with a gold nuclei. 



Projectile 



Energy 



"Inverse kinematic reaction. 
''Fit for generalized formula. 



AA 





0.0 


7.2 ± 1.1 


11+'- 


0.1 


8.7 ± 1.0 


n+'- 


0.18 


10.4 ± 1.2 


7r+'- 


0.3 


12.2 ± 1.3 


P 


0.2 


8.6 ±1.1 


P 


0.49 


14.9 ± 1.8 


P 


0.8" 


17.0 ± 1.4 


P 


0.76 


17.8 ±2.3 


P 


0.8 


18.0 ±2.3 


P 


1.0 


20.7 ±2.8 


P 


1.2 


23.2 ±3.2 


P 


1.6 


26.5 ±3.4 


P 


2.6' 


30.1 ± 1.9 


P 


3.0 


30.6 ±4.3 


P 


11.5 


30.3 ±4.2 


P 


800.0 


26.1 ±3.7 


P 


0.0' 


17.3 ±1.3 


P 


0.0 


20.0 ±0.8 



Ref. 



I|3| 

M 

this work 



distribution of products of gold projectile fragmentation 
on C and Al targets To compare with our results 

we present Sstab obtained for the inverse kinematic Au + 
p reaction at 800 MeV |2^, ^ and for energetic protons 
1^, Besides this, we have also plotted curves derived 
for the Au target from a general formulae describing the 
Zp{A) path for products of various medium and heavy 
target fragmentations induced with protons at the frag- 
mentation limit and with stopped antiprotons [ p7| . 
The lower part of Fig. ^ shows the dependence of the 
charge distribution width on the product mass. 

As can be seen from Fig. |^, the relative course of Zp 
and Zp for stopped antiprotons changes now much more 
dramatically than in Fig. ^. The Sstab calculated for 
more general Y{A, Z) models ^ is smooth due to 
the broad A range fitting. Zp fitted for 11.5 and 300 
GeV protons lies very close to this curves in the 
evaporation residues region. Surprisingly, Zp derived for 
the inverse reaction at smaller energy extends farther to- 
wards the neutron-deficient nuclei for heavy products. In 
the fission region the situation is reversed except for the 
lightest products. The curve plotted for 3 and 6 GeV 
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protons lies closer to the valley of stability, for 1 GeV 
protons only the fission region is represented as there are 
no fit parameters given in Ref. [ p7| . 

When Y{A, Z) models arc applied to shorter mass 
ranges a rather non-continuous behavior with segments 
of rapidly changing position and orientation is observed. 
This happens both for protons and heavy ions 

pO| , 5^ investigated with 7-ray spectroscopy tech- 
nique and for heavy ion reactions on Au studied with the 
inverse kinematic technique p9[ . Obviously, this situa- 
tion cannot be ascribed merely to the uncertainties of the 
experimental data, even in the worst cases. In our case, 
the experimental data distribution in the N — Z plane 
shown in Fig. Q strongly favors the segmentation of the 
evaporation region in fitting. 

Generally, lower energy reactions lead to Zp running 
closer the valley of stability, with passages to the neutron- 
rich side for fission fragments. Antiproton data show a 
peculiar tendency: although Zp lies quite away from the 
valley of stability for evaporation residues, it does not 
reach such a neutron-deficient region as the energetic pro- 
tons. Such a behavior may be partially explained by the 
influence of shell effects observed in antiproton distribu- 
tion for N/Z crossings at 106/76, 82/64 and 50/40. The 
Y{A, Z) yield reaches local maxima in these regions, the 
most probable Z goes towards the more neutron-deficient 
nuclei and the charge dispersion becomes broader, as can 
be seen from Fig. 0. 

The width of the Z distribution, cr, was found to de- 
crease smoothly with decreasing A when a generalized 
approach is used for protons or antiprotons 
On the contrary, results of fitting within shorter mass 
regions are again inconsistent with compilations using 
broad mass regions, with quite small widths for evapo- 
ration residues and with a large scatter of a for fission 
fragments 
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. For antiprotons stopped in Au the 
charge width is rather small in the evaporation region, es- 
pecially in comparison with 800 MeV/nucleon Au on H 
data. On the other hand, products of fission induced by 
stopped antiprotons are distributed quite broadly, simi- 
larly to the proton reaction products. 

From a methodological point of view, results on Zp [A) 
and (J [A) fitted in different ways are not consistent, even 
for protons at similar energies. The Y{A, Z) modeling 
applied for wide mass regions may be reasonable for lim- 
ited experimental data and for generalization purposes, 
however, this approach washes out any possible feature of 
more discrete nature. Hence, the division of input data to 
some A subregions should work better in detailed studies, 
especially for lower excitation reactions. 



charge yield models were fitted only for a part of these 
targets, the parameters of the mass yield distributions for 
the rest were evaluated on the basis of summed direct ex- 
perimental Y{A) yields. However, either for the former 
or the latter results, the yields for the heaviest nuclei, 
close to the target, are underestimated since a significant 
part of the total Y(A) is hidden in non-detectable iso- 
topes. To take this effect into account, we recalculated 
AA values obtained for other targets in the way as it 
was done for Au (see Sec. VA). The results are listed 
in Tabic IV. The removed mass increases with increas- 
ing target mass, as illustrated in Fig. ||. Such behavior 
is consistent with the simple geometrical picture of an 
excitation energy proportional to the number of partic- 
ipating nucleons ||^, hence to the volume of the nucleus 
bombarded with annihilation mesons. 



TABLE IV: Characteristics of mass yield distribution after 
antiproton stopping in different targets. 



Target 



Ref. 



««Ru 
9«Zr 

Cd 

nat^g 

"OTe 
""'Ba 

144 r 



148 J 

i^^Ho 

176Yb 



13.9 ± 1.8 

15.8 ± 1.9 
16.2 ±2.0 
16.2 ± 2.1 

16.4 ±2.3 
17.7 ± 1.8 

18.9 ± 1.6 
17.9 ±2.1 

17.5 ± 1.9 
19.4 ±2.0 
19.2 ± 2.2 
21.7 ±2.5 
21.1 ± 2.0 
22.4 ±2.2 
20.0 ±0.8 



Using the AA value obtained for the Au target we 
may estimate the mean thermal excitation energy of the 
decaying system. The compilation of the measured par- 
ticle emission ^ || gives 5.4 nucleons ejected in the 
cascade±preequilibrium stages through n, p, d, t, '^He 
and ^He ejectiles. Hence we have on the average 14.6 
evaporated nucleons and assuming 8 MeV separation en- 
ergy and 3 MeV kinetic energy per nucleon § leads to 
(161 ±23) MeV stored in the thermalized system. Such a 
result compares nicely with the value of (183 ± 21) MeV 
derived from the measurements of the spectra of neutrons 
and of light charged particles S . 



B. Antiprotons stopped in various targets 

There were many other targets irradiated with low en- 
ergy antiprotons from LEAR |l2|, |l|, |l|, |l5[ |l|, 
review of some results of these experiments wil 
us to look closer at antiproton induced reactions. 
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Mass- 



C. Odd-even effects 

Data on odd-even and shell effects observed in the 
y(A, Z) yield distribution are rather rarely discussed in 
the yield modeling context |21, 37 1. Their influence on 
the yields is difficult to observe if the experimental data 
set is limited and errors are of the order of the possible 
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FIG. 8: Average mass removed from different targets after 
stopped antiproton absorption. Results from earlier experi- 
ments studying mass yield curves ^ ^ ^[ ^| ^] , are 
recalculated after yield correction for the heaviest products. 



odd-even correction. The conclusion of Rudstam po[ , 
looking for a general formula predicting cross sections 
for p and a induced reactions, was that there is no need 
to introduce such a correction as the experiment/model 
yield ratios for various N and Z combinations do not show 
any clear correlation with the nucleon number evenness. 
Later, Silberberg and Tsao found a moderate effect, 
modeled with factors equal to 1.25, 0.9, 1.0 and 0.85 for 
even-even, odd-N, odd-Z and odd-odd (N,Z) pairs, re- 
spectively. 

Since off-line nuclear spectroscopy was applied to study 
the Y{A, Z) distribution the odd-even effect was observed 
in reactions induced with stopped antiprotons jl^, p^ . 
Results for lighter targets (^'^•^^•^^Mo, "^*Ba) are consis- 
tent, with an 18-26 % correction for odd-Z nuclei and 
a 32-34 % correction for odd-N nuclei (and the sum of 
these values in the odd-odd case). Corresponding values 
fitted for ^®^Ho are not so evident, the yield of odd-N 
nuclei is strongly reduced (by w 66%) whereas there is no 
need to correct the odd-Z results. Because these fits were 
made simultaneously for the whole heavy residue region, 
no dependence on the emitted number of nucleons (hence 
excitation) was studied, also no indication for any shell 
effects was reported. 

Using the heavy and fissionable gold nuclei to absorb 
antiprotons we have the opportunity to investigate the 
odd-even and shell effects in a wide evaporation and fis- 
sion products mass range. In Fig. ^ we present three 
fitted yield distributions, as function of mass, charge and 
neutron number of the products. Corrections fitted for 
odd-N and odd-Z yield were shown in Table || with pa- 
rameters ai4 and ois, respectively. They seem to be 
largest for the medium AA region (161 < A < 182), 
since between the heavier residues only the odd-N Pt 
isotopes exhibit a clear yield reduction of about 30%. 
A small odd-N effect (< 10%) seems also to appear for 
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FIG. 9: Yield of residues after gold fragmentation by stopped 
antiprotons: (a) as a function of the atomic mass A, (b) as 
a function of the neutron number N, (c) as a function of the 
atomic number Z. 



the other heavy even-Z products. Lighter evaporation 
residues are produced more uniformly over changing nu- 
cleon numbers evenness, although for products close to 
the closed iV = 82 shell (143 < ^ < 149) the even-N 
isotopes are strongly favored. We have not been able to 
study odd-even effects for the lightest evaporation and all 
fission products because of their small cross sections with 
large relative uncertainties and since the experimental 
data are there much less numerous. It should be stated 
that the correction for odd-odd nuclei used here slightly 
differs from those applied before for stopped antiprotons 
p^ , as we use the multiplicative form (Eq. (||)) in- 
stead of a correction factor equal to 1 — P„ — Pp, where 
Pn and Pp are parameters fitted for odd-N and odd-Z 
nuclei, respectively. Since we observe that the odd-even 
effect is stronger when both N and Z are odd or even, 
the multiplication is more adequate. 

The use of odd-even corrections for stopped antiproton 
reactions strongly improves the fit, with reduction by 
more than one order of magnitude. Thus it cannot be 
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treated as a trivial improvement via adding another pa- 
rameter. Moreover, since corrections obtained for differ- 
ent mass regions are consistent, therefore, taking into ac- 
count previous resuhs obtained for hghter targets |jl^, |l^ 
the inclusion of such a component is unavoidable in cor- 
rect modeling of the data coming from experiments with 
stopped antiprotons. Its strength, more pronounced than 
in the corresponding energetic proton data at 800 MeV 
p9| , |30t , is one of the most distinct features of the an- 
tiproton absorption induced reactions. 

VI. CONCLUSIONS 

The independent and/or cumulative yields for 114 iso- 
topes produced after absorption of stopped antiprotons 
in gold were measured by using the off-line 7-ray spec- 
troscopy technique. On this basis, with the help of a phe- 
nomenological model, the whole yield distribution was 
extracted for residues ranging from the target mass mi- 
nus one down to the light fission products with mass > 20. 
The fission probability was estimated to be (3.8 ± 0.5)%, 
in agreement with the results of measurements using on- 
line techniques. 

An average thermal excitation energy, gained by the 
Au nucleus after p annihilation, was shown to be quite 
similar to that of 1 GeV protons, although the fission 
probability for such protons is almost twice as large. 
Moreover, the inspection of the yield distribution over 



the A-Z plane indicates a fairly peculiar character of the 
reaction induced with low energy antiprotons. The most 
probable Z course is quite different, lying closer to the 
stability line and exhibiting a more complex shape. Fur- 
thermore, the charge dispersion over Z does not compare 
with that observed for 0.8 GeV protons, being almost 
twice as narrow. 

The average mass removal observed for various targets 
reacting with stopped antiprotons rises linearly with in- 
creasing target mass. This behavior is consistent with in- 
beam studies of the light particle emission. derived 
from mass yield data helps to complete such measure- 
ments, unable to detect charged particles of the lowest 
energy. 

A clear odd-even and some shell effects distinguish ev- 
idently the reaction with stopped antiproton from those 
induced with energetic protons. For the first time the 
dependence of this phenomenon on the residue mass was 
studied. The strength of such effects seems to diminish 
with the excitation energy, although for long evaporation 
chains and fission products it may be unobserved due to 
scarce and uncertain data. 
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